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Hold-in, pull-in, and lock-in ranges of PLL circuits: 
rigorous mathematical definitions 
and limitations of classical theory. 

Leonov G.A., Kuznetsov N.V., Yuldashev M.V., Yuldashev R.V. 


Abstract —The terms hold-in, pull-in (capture), and 
lock-in ranges are widely used by engineers for the 
concepts of frequency deviation ranges within which 
PLL-based circuits can achieve lock under various ad¬ 
ditional conditions. Usually only non-strict definitions 
are given for these concepts in engineering literature. 
After many years of their usage, F. Gardner in the 2nd 
edition of his well-known work, Phaselock Techniques , 
wrote u There is no natural way to define exactly any unique 
lock-in frequency” and “despite its vague reality, lock-in range 
is a useful concept Recently these observations have 
led to the following advice given in a handbook on 
synchronization and communications “ We recommend 
that you check these definitions carefully before using them ” 11 , 
p.49]. In this survey an attempt is made to discuss and 
fill some of the gaps identified between mathematical 
control theory, the theory of dynamical systems and 
the engineering practice of phase-locked loops. It is 
shown that, from a mathematical point of view, in 
some cases the hold-in and pull-in “ranges” may not 
be the intervals of values but a union of intervals and 
thus their widely used definitions require clarification. 
Rigorous mathematical definitions for the hold-in, pull- 
in, and lock-in ranges are given. An effective solution 
for the problem on the unique definition of the lock-in 
frequency, posed by Gardner, is suggested. 

Index Terms —Phase-locked loop, nonlinear analysis, 
analog PLL, high-order filter, local stability, global 
stability, stability in the large, cycle slipping, hold-in 
range, pull-in range, capture range, lock-in range, defi¬ 
nition, Gardner’s problem on unique lock-in frequency, 
Gardner’s paradox on lock-in range. 


I. Introduction 

T HE phase-locked loop based circuits (PLL) are 
widely used in various applications. A PLL is es¬ 
sentially a nonlinear control system and its nonlinear 
analysis is a challenging task. Much engineering writing 
is devoted to the study of PLL-based circuits and the 
various characteristics for their stability (see, e.g. a rather 
comprehensive bibliography of pioneering works in !)• 
An important engineering characteristic of PLL is a set of 
parameters’ values for which a PLL achieves lock. In the 
classical books on PLLs i published in 1966, such 
concepts as hold-in, pull-in, lock-in, and other frequency 


ranges for which PLL can achieve lock, were introduced. 
They are widely used nowadays (see, e.g. contemporary 
engineering literature M and other publications). Usu¬ 
ally in engineering literature only non-strict definitions 
are given for these concepts. F. Gardner in 197^] in the 
2nd edition of his well-known work, Phaselock Techniques , 
formulated the following problem |[9, p.70] (see also the 3rd 
edition [6] p. 187-188]): “ There is no natural way to define 
exactly any unique lock-in frequency ”. The lack of rigorous 
explanations led to the paradox: “ despite its vague reality, 
lock-in range is a useful concept ” |9, p.70]. Many years of 
using definitions based on the above concepts has led to 
the advice given in a handbook on synchronization and 
communications, namely to check the definitions carefully 
before using them 1, p.49]. 

In this paper it is shown that, from a mathematical 
point of view, in some cases the hold-in and pull-in 
“ranges” may be not intervals of values but a union of 
intervals, and thus their widely used definitions require 
clarification. Next, rigorous mathematical definitions for 
the hold-in, pull-in, and lock-in ranges are given. In ad¬ 
dition we suggest an effective solution for the problem of 
the unique definition of the lock-in frequency, posed by 
Gardner. 


II. Classical nonlinear mathematical models of 
PLL-based circuits in a signal’s phase space 


In classical engineering publications various analog 
PLL-based circuits are represented in a signal’s phase 
space (also named frequency-domain 10 p.338]) by the 
block diagram shown in Fig. [l] 
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Fig. 1. PLL-based circuit in a signal’s phase space. 

1 A year later, in 1980, F.Gardner was elected IEEE Fellow “for 
contributions to the understanding and applications of phase lock 
loops”. 
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Considering the corresponding mathematical model: the 
Phase Detector (PD) is a nonlinear block; the phases 
^ 1,2 (f) of the input (reference) and VCO signals are PD 
block inputs and the output is a function c p(0 A (t )) = 
ip(Qi(t) — # 2 (t)) named a phase detector characteristic , 
where 

0 A (t) = Oi{t) ~ (t), (1) 

named the phase error. The relationship between the input 
(p{0 A (t)) and the output g(t) of the linear filter (Loop 
filter) is as follows: 

x = Ax + bcp(0 A (t)) : g(t ) = c*x + hip(0 A (t)), (2) 

where A is a constant matrix, x(t) G M n the filter state, 
x(0) the initial state of filter, b and c constant vectors, and 
h a number. The filter transfer function has the form0 

H(s) = —c*(A - siy'b + h. (3) 

A lead-lag filter [t] (usually 0) = —c*A~ 1 b + h = 1, 
but H{ 0) can also be any nonzero value when an active 
lead-lag filter is used), or a PI filter (H( 0) is infinite) is 
usually used as the filter. The solution of 0 with initial 
data x(0) (the filter output for the initial state x(0)) is as 
follows: 

g(t,x 0 ) = a 0 (t,x(0)) + f j(t- t)<p(0 a (t))(1t + h<p(0 A (t)), 
0 

(4) 

where j(t — r) = c*e A( d~ r ^b is the impulse response 
function of the filter and <ao(£,x(0)) = c*e At x( 0) the zero 
input response (natural response, i.e. when the input of 
the filter is zero). The control signal g(t) adjusts the VCO 
frequency to the frequency of the input signal: 

0 2 (t) = U} 2 (t) = t4 ree + Lg(t), (5) 

where c<4 ree is the VCO free-running frequency (i.e. for 
g{t) = 0) and L the VCO gain. Nonlinear VCO models can 
be similarly considered, see, e.g. [l2 , [l3]. The frequency 
of the input signal (reference frequency) is usually assumed 
to be constant: 


0i{t) = cji(t) = cui. (6) 


The difference between the reference frequency and the 
VCO free-running frequency is denoted as c<;^ ee : 


c4 ee = cji - 4 ree . 


(7) 


By combining equations Q, Q, and 0-0 a nonlinear 
mathematical model in the signal’s phase space is obtained 
(i.e. in the state space: the filter’s state x and the difference 
between the signal’s phases 0 A )\ 


x = Ax + b(p(0 A ), 

6 A = a4 ee — Lc*x — Lh<p(Q A ). 


Nowadays nonlinear model (|8| is widely used (see, e.g. 7 , 
14 , [l5]) to study acquisition processes of various circuits. 


The model can be obtained from the corresponding model 


2 In the control theory t he t ransfer function is often defined with 
the opposite sign (see, e.g. 11 ): H(s) = c*(A — sl)~ 1 b — h. 


in the signal space (called also time-domain 
by averaging under certain conditions 11 


10 


16 


p.329]) 

|p], a 

rigorous consideration of which is often omitted (see, e.g. 
classical books [4j p.12,15-17], [3, p.7]) while their violation 


may lead to unreliable results (see, e.g. 20 , |21|). 


Usually the PD characteristic is an odd function (e.g. 
a PD realization such as a multiplier, JK-flipflop, EXOR, 
PFD, and other elements 0)- Note that the PD charac¬ 
teristic <p(0 A ) depends on the waveforms of the considered 
signals [l8], 1 19 ). For the classical PLL with sinusoidal sig¬ 
nals and a two-phase PLL we have ip{0 A ) = \ sin (0 A ), for 
the classical BPSK Costas loop with ideal low-pass filters 
and a two-phase Costas loop we have ip(0 A ) = | sin(2#A)- 
Classical PD characteristics are bounded piecewise 
smooth 27r periodic function^] 

<p{0 A (t) + 27rfc) = <p(0 A (t)), Vfc = 0,1, 2... 


Thus, it is convenient to assume that 0 A mod 27r is a 
cyclic variable, and the analysis is restricted to the range 
Of #a(0) G [—7T,7r). 

For the case of an odd PD characteristic]^] system 0 is 
not changed by the transformation 

(c 0 % ee ,x(t),e A (t)) ^(-j™,-x(t),-e A (t)). (9) 

Property 0 allows the analysis of system (]8| with only 
^ q an q introduces the concept of frequency deviation 

I^A e | = l^ 1 -^2 ee |- 


III. Locked state 

The locked states (also called steady states) of the 
model in the signal’s phase space must satisfy the following 
conditions: 

• the phase error 6 A is constant, the frequency error 0 A 
is zero; 

• the model in a locked state approaches the same 
locked state after small perturbations (of the VCO 
phase, input signal phase, and filter state). 

The locally asymptotically stable equilibrium (stationary) 
points of model (|8|: 

0 A (t) = 0 eq + 27T&, x{t) = X eq , (10) 

are locked states, i.e. satisfy the above condition^] 

Considering the case of a nonsingular matrix A (i.e. the 
transfer function of the filter does not have zero poles), 
the equilibria of (|8| (stationary points) are given by the 
equations 

, .free . .free 

, n /n \ _ 

^ eq> L(-cM- 1 6 + /i) LH{ 0)’ 

free ' 

X eq = -A^b^Oeq) = 


3 If p(0 A {t)) has another period (e.g. 7r for the Costas loop 
models), it has to be considered in the further discussion instead 
of 27r. 

4 There are examples of non odd PD characteristics, where (19} doe s 
not hold true (see, e.g. BPSK Costas loop with sawtooth signals [18] , 

19] and others). 

0 It can be proved that if the filter is controllable and observable, 
then only equilibria satisfy locked state co ndit ions, i.e. the filter state 
x{t) must be constant in the locked state 11 . 
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Thus, the equilibria can be considered as a multiple-valued 
function of variable co^ ee : [x eq {i <;^ ee ), 6 eq (uj% ee )) . From the 
boundedness of the PD characteristic c p(6 eq ) it follows that 
there are no equilibria for sufficiently large IcJa 66 ! (see 
Figs. | and §. 


IV. Engineering definitions of stability ranges 


The widely used engineering assumption (see Viterbi’s 
pioneering writing [4j p. 15]) is that the zero input response 
of filter <To(t,xo) does not affect the synchronization of 
the loop. This assumption allows the filter state x(t ) 
to be excluded from the consideration and a simplified 
mathematical model of PLL-based circuit in the signal’s 
phase space to be obtained from Q and © ( see > e.g. 
p.17, eq.2.20] for h = 0 and [3j p.41, eq.4-26] for 7 = 0): 

6 A =c o% ee - Lh(t - t)lp(6a(t))&t - Lh(p(6 A {t))- (12) 
Jo 


For an example of this one-dimensional integro-differential 
equation the following intervals ( @) are defined: the 

hold-in range includes |co^ ee | such that model © has 
an equilibrium 6 A {t) = 0 eq , which is locally stable (local 
stability, i.e. for some initial phase error 0 A { 0)); the pull-in 
range includes \u% ee \ such that any solution of model © 
is attracted to one of the equilibria 0 eq (global stability, i.e. 
for any initial phase error 0 a (0)). Thus, the block diagram 
of the loop in Fig. [l] is usually considered without initial 
data x(0) and 0 a (0) (see, e.g. @ p- 17 > fi g- 2 - 3 ])- 

Viterbi [ 4 ] explains the above assumption for the stable 
matrix A, but considers also various filters with marginally 
stable matrixes (e.g. a filter - perfect integrator, where 
A = 0). At the same time, even for a stable matrix A, the 
initial filter state x(0) and ao(t,xo) may affect the acqui¬ 
sition process and stability ranges (see, e.g. corresponding 
examples for the classical PLL 20 and Costas loops 

m 


21 


While the above assumption allows introduction of 
the above one-dimensional stability sets, defined only by 
|^ ee |, for rigorous study the multi-dimensional stability 
domains have to considered, taking into account x(0), and 
their relationships with the classical engineering ranges 
have to be explained. In (6j p. 187] it is noted that the 
consideration of all state variables is of utmost importance 
in the study of cycle slips and the lock-in concept. 


V. Rigorous definitions of stability sets 

The rigorous mathematical definitions of the hold-in, 
pull-in, and lock-in sets are now given for the nonlinear 
mathematical model of PLL-based circuits in the signal’s 
phase space (J8| and corresponding nontrivial examples are 
considered. 


A. Local stability and hold-in set 

We now consider the linear iz at iorj^] of system (pi along 
an equilibrium (x eq ,6 eq ). Taking into account plj ) and 
:= d(p(0)/d6, the linearized system is as follows: 

( ^ ^ ( A bip (6 eq ) \ f x — x eq \ , s 

\o A J V - Lc * - ) \o A -o eq ) [i6) 

The characteristic polynomial of linear system (pi) can be 
written (using the Schur complement, e.g. llj) in the 
following form: x( 5 ) = ( — Lh(p'(0 eq ) — s + Lc*(A — 
8/)- 1 6(^ / (0 eg ))det(A-5/), or can be expressed in terms of 
the filter’s transfer function H(s) = ^y, where a(s) and 
d(s) are polynomials: 

x(s) = -(sd(s) + a(s)L<p'(0 e q)). (14) 


The characteristic polynomial corresponds to the denom¬ 
inator of the closed loop transfer functiorQ 

To study the local stability of equilibria ( pTj ), it is nec¬ 
essary to check whether all the roots of the characteristic 
polynomial (14) for the linearization of model © along 
the equilibria (i.e. the poles of the closed loop transfer 
function) have negative real parts. For this purpose, at 
the stage of pre-design analysis when all parameters of the 
loop can be chosen precisely, the Routh-Hurwitz criterion 


and its analogs (see, e.g. Kharitonov’s generalization 30 


for interval polynomials) can be effectively applied. At 
the stage of post-design analysis when only the input and 
VCO output are considered and the parameters are known 
only approximately, various frequency characteristics of 
the loop (see, e.g. Nyquist and Bode plots) and the 
continuation principle can be used (see, e.g, ©. HP- 

If the PD characteristic is an odd function and hence 
( p'(0 eq ) is an even function, from © we conclude that 

1) there are symmetric equilibria: 

p eg (u£nM^A e )) = (-x eg (-^ ee ),-0 e9 (-^ ee )), 

2) these symmetric equilibria are simultaneously stable 
or unstable. 

The same holds true for nonstationary trajectories. 


Definition 1. A set of all frequency deviations |cj^ ee | such 
that the mathematical model of the loop in the signal’s 
phase space has a locally asymptotically stable equilibrium 
is called a hold-in set Lipoid-in- 


Thus, a value of frequency deviation belongs to the hold- 
in set if the loop re-achieves locked state after small per¬ 
turbations of the filter’s state, the phases and frequencies 
of VCO and the input signals. This effect is also called 
steady-state stability. In addition, for a frequency deviation 
within the hold-in set, the loop in a locked state tracks 


6 Here it is assumed that the PD characteristic is smooth 

at the point 6/± = 6 eq . However, there are PLL-based circuits 
with nonsmooth or discontinuous PD characteristics (see, e.g. the 
sawtooth PD characteristic for PLL IS] , the model of QPSK Costas 
loop [25], and some others [26 -[28]). In such a case care has to be 
taken of the definition of solutions, the linearizatio n of the model and 
the analysis of possible sliding solutions (see, e.g. [29] ). 

7 Consideration of linearized model \13\ allows to avoid the rig¬ 
orous discussion of initial states (x(0), 0/^(0)) related to the Laplace 
transformation and transfer functions 11 . 
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Fig. 2. Phase portrait for u;^ ee from the hold-in range. The 
system’s evolving state over time traces a trajectory (x(t), 0/±(t)). 
Trajectories can not intersect. Unstable equilibrium points, such as 
saddles — black dots, locally asymptotically stable equilibria — 
green dots, any of which has its own basin of attraction (shaded 
domain) bounded by stable saddle separatrices (black trajectories 
going to the saddles). There are initial states and corresponding 
trajectories (see, e.g. dashed trajectory), which are not attracted to 
an equilibrium. 


small changes in input frequency, i.e. achieves a new locked 
state (tracking process ). 

In the literature the following explanations of the hold- 
in range (sometimes also called a lock range 31 p.507], 


32, p.10-2], a synchronization range 33 , a tracking range 
1, p.49]) can be found: “The hold-in range is obtained by 
calculating the frequency where the phase error is at its 
maximum ” [34], p. 171], “The maximum frequency differ¬ 
ence before losing lock of the PLL system is called the hold- 
in range ” [8j p.258]. The following example shows that 
these explanations may not be correct, because for high- 
order filters the hold-in “range” may have holes. 

The following example shows that the hold-in set may 
not include cj^ ee = 0. 


Example 1 (the hold-in set does not contain cj^ ee = 0). 
Consider the classical PLL with the sinusoidal PD char¬ 
acteristic <p(0a) = \ sui(#a), VCO input gain L = 8, and 
the filter transfer function 


H(s) = = 1 + O- 55 

1 ’ d(s) 1 + 0.5s + 0.5s 2 ' 


(15) 


From (11) the following equation for equilibria is obtained: 


1 

2 


sin(# e< j) 



(16) 


Applying the Routh-Hurwitz stability criterior f] t o the 
denominator of the closed loop transfer function (14) 


s 3 + s 2 + s( 2 + 4 cos (0 eq )) + 8 cos (0 eq ), 


(17) 


8 For a third-order polynomial x( s ) — &3-s 3 + « 2 <s 2 + a\s + ao, 
all the roots have negative real parts and the corresponding linear 
system is asymptotically stable if < 21 , 2,3 > 0 and 0 , 20-1 > 0300 . For 
X(s) = 04 s 4 T 03 s 3 + 02 S 2 + ois + ao, all the coefficients must satisfy 
01 , 2 , 3,4 > 0, and 0302 > 0401 and 030201 > 040 2 + o|oq. 



Fig. 3. Phase portrait for o;^ ee outside the hold-in range: there are 
no locally stable equilibria. 


the following conditions are obtained: 

cos(0 eq ) > 0, (2 + 4cos(6» eg )) > 0, 

(2 + 4cos(# e9 )) > 8cos(# e< j). 

Then 0 < cos(9 eq ) < \. and for the locked state the steady- 
state phase error (i.e. corresponding to an equilibrium) is 
obtained 

^e 9 e(--,--) u (-,-). (19) 

Taking into account (Jl6| ; ( [l9| ) ; one obtains the hold-in set 

|^ ee |G(2y3,4). (20) 


The next example shows that the hold-in set may not 
actually be a range (i.e., an interval) but a union of 
intervals, one of which may include Go^ ee = 0. 


Example 2 (the hold-in set is a union of disjoint intervals, 
one of which contains co^ ee = 0). Consider the classi¬ 
cal PLL with the sinusoidal PD characteristic </?(#a) = 
4sui(#a); the VCO input gain L = 80, and the filter 
transfer function 


1 + 0.25s + 0.5s 2 
1 + 2s + 2s 2 + 2s 3 ' 


( 21 ) 


From ( [IT] ) the following equation for the equilibria is 
obtained: 

2 sin (M = ( 22 ) 


An equilibrium is asymptotically stable if and only if all the 
roots of polynomial (14); 


5(1 + 25 + 25 2 + 25 s ) + K( 1 + 0.255 + 0.55 2 ) = 

25 4 + 25 s + 5 2 (2 + 0.5 K) + 5(1 + 0.25i^) + K, (23) 
K = Lip\6 eq ) = 40 cos (0 eq ) 


have negative real parts. Using the Routh-Hurwitz crite¬ 
rion, we obtain 


2 + 0.5 K > 0, 1 + 0.25 K >0, K > 0, 

2(2 + 0.5 K) > 2(1 + 0.25 K), (24) 

2(2 + Q.5K)(1 + 0.25i^) > 2(1 + 0.25if) 2 + 2 2 ^- 
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From these inequalities we have 

K = 40 cos(9 eq ) e (0,12 - 8\/2) U (12 + 8\/2, oo), 

0 e q e (-G —1.5536) U (-0.9486,0.9486) U (1.5536, (-). 

^ 4 

(25) 

Note that for other values of 0 eq at least one root of 
the polynomial (23) has a positive real part, making the 
corresponding equilibrium unstable. Combining (22) and 
(25) ; we obtain the hold-in set 


\J£ ee \ € [0,32.5) U (39.9942,40). 


(26) 


Note that in this case, for the values of the VCO input gain 
L > 24 + 16\/2 the hold-in set is always a union of disjoint 
intervals. For L = 80 the simulation results of transition 
process in Simulink mode^\in Fig.are shown in Figs. M 
for the initial data (x(0) = (0; 0; 0.9990), #a(0) = 1.5585) 
and various uj^ e . 



free-running 

frequency 


Fig. 4. MatLab Simulink: the signal’s phase space model of the 
classical PLL 


Related discussion on the frequency responses of loop 
with high-order filters can be found in [6j p.34-38, 52-56]. 

Remark 1. For the first order filters, the set f ^hold-in I s 
an interval |cj^ ee | < ujh- For higher order filters, the set 
hold-in may be more complex. Thus, from an engineering 
point of view, it is reasonable to require that = 0 

belongs to the hold-in set and to define a hold-in range as 
the largest interval [0,^) from the hold-in set 

[ 0 , idh) ^ ^hold -in 

such that a certain stable equilibrium varies continuously 
when is changed within the rang^J Here Uh is called 
a hold-in frequency (see @ P-38J). 

9 Following the above classical consideration, the filter is often 
represented in MatLab Simulink as the block Transfer Fen with zero 
initial state (see, e.g. [35] — [39] ). It is also related to the fact that the 
transfer function (from (p to g) of linear system pi is defined by the 
Laplace transformation for zero initial data x(0j1 = 0. In Fig. [4] we 
use the block Transfer Fen (with initial states) to take into account 
the initial filter state :r(0); the initial phase error $a( 0) can be taken 
into account by the property initial data of the Intergator blocks. 
Note that the corresponding initial states in SPICE (e.g. capaci tor’ s 
initial charge) are zero by default but can be changed manually 40 . 

10 In general (when the stable equilibria coexist and some of them 
may appear or disappear), the stable equilibria can be considered as a 
multiple-valued function of variable o;^ ee , in which case the existence 
of its continuous singlevalue branch for |a;^ ee | € [0,^) is required. 


Remark 2. In the general case when there is no symmetry 
with respect to c (see the hold-in set need not be 
symmetric and the set uo^ ee G Tthoid-in must be considered 
in Definition [7} 



Fig. 5. o;^ ee = 3: stable locked state exists. 




Fig. 7. ca^ ee = 39.997: stable locked state exists. 


B. Global stability (stability in the large) and pull-in set 

Assume that the loop power supply is initially switched 
off and then at t = 0 the power is switched on, and assume 
that the initial frequency difference is sufficiently large. 
The loop may not lock within one beat note, but the 
VCO frequency will be slowly tuned toward the reference 
frequency (acquisition process ). This effect is also called a 
transient stability. The pull-in range is used to name such 
frequency deviations that make the acquisition process 
possible (see, e.g. explanations in |3j p.40], 0 P-61])- 
To define a pull-in range (called also a capture range 
41 , an acquisition range (33l p.253]) rigorously, consider 


first an important definition from stability theory. 

Definition 2. If for a certain uj^ e any trajectory of sys¬ 
tem <® tends to an equilibrium, then the system with such 
is called globally asymptotically stable (see Fig. [#|). 
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Fig. 8. Phase portrait for cu^ ee from the pull-in range: any trajectory 
is attracted to an equilibrium (equilibria: green—stable and black- 
unstable circles); for a sufficiently large initial state of the filter, cycle 
slipping is possible (see, e.g. dashed trajectory). 


We now consider a possible rigorous definition. 

Definition 3. A set of all frequency deviations |co^ ee | such 
that the mathematical model of the loop in the signal’s 
phase space is globally asymptotically stable is called a pull- 
in set Qpuii -in • 

Remark 3. In the general case when there is no symmetry 
with respect to the set uj^ e G Q pu ii-in has to be 

considered in Definition [3j 

Remark 4. The pull-in set is a subset of the hold-in set: 
^pull-in C Dhoid-in, and need not be an interval. From an 
engineering point of view, it is reasonable to require that 
^free _ q ]j e i on g S pull-in set and to define a pull-in 

range as the largest interval [0,cj p ) from the pull-in set: 

[0,CiJp) C Hpull-im 

where uj v is called a pull-in frequency (see § p-40]). 
Remark 5. If all possible states of the filter are bounded: 


x G X rea i (e.g. X rea i — {x : c m i n < \x\ < c max )f 


by the design of the circuit (e.g. capacitors have limited 
maximum and minimum charges, the VCO frequency is 
limited etc.), then in the definition of pull-in set it is 
reasonable to require that only solutions with x(0) G X rea i 
tend to the stationary set. Trajectories, with initial data 
outside of the domain defined by x(0) G X rea i (here the 
initial phase error 0 a (0 ) can take any value), need not tend 
to the stationary set. 


For the model without filter (i.e. H(s) = const ) the 
pull-in set coincides with the hold-in set. The pull-in 
set of PLL-based circuits with first-order filters can be 
estimated using phase plane analysis methods [42], 43 , 
but in general its rigorous study is a challenging task [4 , 


12 


17 


44 


45 


For the case of the passive lead-lag filter H(s) = 
i+s(ti+ T 2 ) ? a recent work 12, p. 123] notes that “the de¬ 
termination of the width of the capture range together 
with the interpretation of the capture effect in the second 
order type-I loops have always been an attractive theoretical 


problem. This problem has not yet been provided with a 
satisfactory solution ”. At the same time in 11 , 46-48 


it is shown that the basin of attraction of the stationary set 
may be bounded (e.g. by a semistable periodic trajectory, 
which may appear as the result of collision of unstable and 
stable periodic solutions), and corresponding analytical 
estimations and bifurcation diagram are given. 

Note that in this case a numerical simulation may give 
wrong estimates and should be used very carefully. For 
example, in [40] the SIMetrics SPICE model for a two- 
phase PLL with a lead-lag filter gives two essentially 
different results of simulation with default “auto” sampling 
step (acquires lock) and minimum sampling step set to Ira 
(does not acquire lock — such behaviour agrees with the 
theoretical analysis). The same problems are also observed 
in MatLab Simulink 20 , 21 , 49 , see, e.g. Fig. [9] These 
examples demonstrate the difficulties of numerical search 
of so-called hidden oscillations 48 , 50 , 51 , whose basin 


of attraction does not overlap with the neighborhood 
of an equilibrium point, and thus may be difficult to 
find numericall^P') In this case the observation of one or 
another stable solution may depend on the initial data and 
integration step. 



Fig. 9. Simulation of two-phase PLL described by Fig. ^]or model 
|Hn = 0.0448, T2 = 0.0185, A = 6=1- 

c = h = n?= §siii(0 A ); wi = 10000, u> f 2 ree = 

10000 — 178.9, L = 500. Filter output g(t) for the initial data 
xo = 0.1318, #a(0) = 0 obtained for default “auto” relative tolerance 
(red) — acquires lock, relative tolerance set to “le-3” (green) — does 
not acquire lock. 


S. Goldman, who has worked at Texas Instruments over 
20 years, notes that PLLs are used as pipe cleaners for 
breaking simulation tools [54] p.XIII]. 

11 In 52] the crash of aircraft YF-22 Boeing in April 1992, caused 
by the difficulties of rigorous analysis and design of nonlinear control 
systems with saturation, is discussed and the conclusion is made that 
since stability in simulations does not imply stability of the phys¬ 
ical control system, stronger theoretical understanding is required 
(see, e.g. similar problem with the simulation of PLL in Fig. [ 9 ]. 
These difficulties in part are related to well-known Aizerman’s and 
Kalman’s conjectures on the global stability of nonlinear control 
systems, which are valid from the standpoint of simplified analy¬ 
sis by the linearization, harmonic balance, and describing function 
methods (note that all these methods are a lso wide ly used to the 
analysis of nonlinear oscillators used in VCO 12], 13]). However the 
counterexamples (multistable high-order nonlinear systems where 
the only equilibrium, which is stable, coexists with a hidden periodic 
oscillation) can be constructed to these conjectures [48| , [53] . 
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While PLL-based circuits are nonlinear control systems Definition HI# 
and for their nonlocal analysis it is essential to apply the 
classical stability criteria, which are developed in control 
theory, however their direct application to analysis of 
the PLL-based models is often impossible, because such 
criteria are usually not adapted for the cylindrical phase 
spacj^j in the tutorial Phase Locked Loops: a Control 
Centric Tutorial 14 , presented at the American Control 
Conference 2002 , it was said that “ The general theory of 


sup |0 a(O) - 0 A (t)\ > 27r, 
t>o 


(29) 


it is then said that cycle slipping has occurred. 

Note that, in general, Definition [4j need not mean 
that finally (after acquisition) condition ( [28] ) can not be 
fulfilled. 

Sometimes, the number of cycle slips is of interest. 


PLLs and ideas on how to make them even more useful 
seems to cross into the controls literature only rarely ”. 

At the same time the corresponding modifications of 
classical stability criteria for the nonlinear analysis of 
control systems in cylindrical phase space were well de¬ 
veloped in the second half of the 20th century, see, e.g. 
29 , 58 - 60 . A comprehensive discussion and the cur¬ 
rent state of the art can be found in 11 . One reason 


why these works have remained almost unnoticed by the 
contemporary engineering community may be that they 
were written in the language of control theory and the 
theory of dynamical systems, and, thus, may not be well 
adapted to the terms and objects used in the engineering 
practice of phase-locked loops. Another possible reason, as 
noted in [61, p.l], is that the nonlinear analysis techniques 
are well beyond the scope of most undergraduate courses 
in communication theory and circuits design. Note that 
for the application of various stability criteria it is often 
necessary to represent system ^ in the Lur’e form: 

A 0 \ { x\ , { b 


0a 


-Lc * 0 / \0 


+ 


-Lh 


£(0a), (27) 


where 
x = x 


x. 


x + A 1 bip(0 eg ), <p(0a) = 7>(#a) - <f(0 eq ), 


\-l 


eq 

(0 eq ) = uj i ^ e L~ 1 (—c*A~ 1 b + h)~ . 

See also discussion of some nonlinear methods for the 


analysis of PLL-based models in recent books 12 

m p . 


13 


C. Cycle slips and lock-in range 

Let us rigorously define cycle slipping in the phase space 
of system 

Definition 4. If 


limsup |0 a(O) — 0 A (t)\ > 27r, 

t —)*-|-oo 


(28) 


Definition 5. If 

2kir < limsu] 

t — Yog 

it is then said that k cycle slips occurred. 


2kir < limsup |$a(0) — 0a(^)| < 2(fc + l)7r, (30) 

t—YOO V 7 


A numerical study of cycle slipping in classical PLL 
can be found in 63 . Analytical tools for estimating the 


number of cycle slips depending on the parameters of the 
loop can be found, e.g. in 


11 , 1581, 64 


The concepts of lock-in frequency and lock-in range 


(called also a lock range 65, p.256], a seize range 66 


p. 138]), were intended to describe the set of frequency 
deviations for which the loop can acquire lock within one 
beat without cycle slipping. In [3, p.40] the following defi¬ 
nition was introduced: “If for some reason, the frequency 
difference between input and VCO is less than the loop 
bandwidth, the loop will lock up almost instantaneously 
without slipping cycles. The maximum frequency difference 
for which this fast acquisition is possible is called the lock-in 
frequency ”. 

However, in general, even for zero frequency deviation 
(uj% ee = 0) and a sufficiently large initial state of filter 
(x(0)), cycle slipping may take place (see, e.g. dashed 
trajectory in Fig. JToJleft). Thus, considering of all state 
variables is of utmost importance for the cycle slip analysis 
and, therefore, the concept lock-in frequency lacks rigor for 


classical simplified model (12) because it does not take into 


account the initial state of the filter. The above definition 
of the lock-in frequency and corresponding definition of 
the lock-in range were subsequently in various engineering 


J68| p.161], 69] p.612], 
T4TP-41, |72l p.24l, 
P-61], [66, p-138], 


73 


75 


it is then said that cycle slipping occurs (see, e.g. dashed 
trajectory in Fig. [#|). 

Here, sometimes, instead of the limit of the difference, 
the maximum of the difference is considered (see, e.g. 44 
P-131]). 

12 For example, in the classical Krasovskii-LaSalle principle on 
global stability the Lyapunov function has to be radially unbounded 
(e.g. V(x,6a) +00 as ||(x,6> A )|| —> + 00 ). While for the application 
of this principle to the analysis of phase synchronization systems 
there are usually used Lyapunov functions periodic in 6/± (e.g. 
V(x, 6 a ) in Remark [8] is bounded for any || (0,0a) 11 T oo), and 
the dis cus sion of this gap is often omitted (see, e.g. patent 15 and 
works |55|-|57l). Rigorous discussion can be found, e.g. in ll]7[29 . 


publications (see, e.g. ||67, p.34-35 
(70j p.532], 0 p.25], 0 p.49], 
p.749], (74) p.56], (54, p.112], 17 
p.576], |J p.258]). 

The loop model © has a subdomain of the phase 
space, where trajectories do not slip cycles (called a lock-in 
domain), for each value of cj^ ee . The lock-in domain is the 
union of local lock-in domains, each of which corresponds 
to one of the equilibria and has its own shape (see, e.g. 
shaded domain in Fig. [lO] left defined by corresponding 
separatrices). The shape of the lock-in domain significantly 
varies depending on co^ ee . In p.50]) a lock-in domain is 
called a frequency lock. Some writers (e.g. 44 p. 132], 76, 
p.355]) use the concept lock-in range to denote a lock-in 
domain. 

In general, taking into account nonuniform behavior of 
the lock-in domain shape, Gardner wrote “There is no 
natural way to define exactly any unique lock-in frequency” 
§ P-70], § p.188]. 
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— (red) negative co A free 

— (black) positive co A free 



Fig. 10. Phase portraits for the classical PLL with the following parameters: H(s) = Ti — 4.48 -10 2 , T 2 = 1.85 -10 2 , 

L = 250, c p(0 /±) = 4 sin(^ z \), and various frequency deviations. Black color is for the system with positive o;^ ee = |u;|. Red is for the system 
with negative a;^ ee = — |u;|. Equilibria (dots), separatrices pass in and out of the saddles, local lock-in domains are shaded (upper black 
horizontal lines is for o;^ ee > 0, lower red vertical lines is for o;^ ee < 0). Left subfig: o;^ ee = 0; middle subfig: cu^ ee = ±65; right subfig: 

U% ee = ±6g> 


Below we demonstrate how to overcome these problems 
and rigorously define a unique lock-in frequency and range. 
We now consider a specific cj^ ee and denote by 


^lock-in a) the corresponding lock-in domain. Such a 
domain exists for any |cj^ ee | G f^hoid-in because at least 
the equilibria are contained in this domain. For a set 
oj% ee G Q we consider the intersection of corresponding 
lock-in domains (see, e.g. the intersections of local lock-in 

domains 


10 


domains for various co^ ee = ±|cj| in Fig. 
shaded both by red vertical and black horizontal lines): 

Dlock-in(ft) = f| Dlock-in^A 66 )- 

Definition 6. A lock-in range is the largest interval [0,cl?/) 
such that for any |cj^ ee | G [0,c^) the mathematical model 
of the loop in the signal’s phase space is globally asymptoti¬ 
cally stable (i.e. [0,u;/) C [0 ,uj p )) and the following domain 

D/oc/c-m(( = ^ lock-ini^'A )* 

\" f r\<u>i 


contains all corresponding equilibria: 

(x eq (^ ee ),0 e g(u} f ^ ee )) e D lock-in((-Ul,UJi)). 

We call such domain Di oc k_ in = Di oc k_i n ((—> ^i)) a uni¬ 
form lock-in domain (uniform with respect to (—a;/,a;/)), 
uoi is called a lock-in frequency (see ! p.40]). 

Various additional requirements may be imposed on the 
shape of the uniform lock-in domain Di oc k_in, e.g. it has to 
contain the line defined by x = 0 (see, e.g. (8j p.258]) or the 
band defined by \x\ < c max . If instead of global stability 
in the definition of the pull-in set we consider stability 
in the domain defined by Weal, then we require that 
the intersection Di oc k_i n ("] Weal contains all corresponding 
equilibria. 


Remark 6. In the general case when there is no symmetry 
with respect to u: A ee we have to consider a unsymmetrical 
interval containing zero in Definition [#| 


Similarly, we can define an extension of the lock-in range: 
^lock-in D [0,cj/), called a lock-in set (however, in general, 
such an extension may be not unique). 

In other words, the definition implies that if the loop 
is in a locked state, then after an abrupt change of u: A ee 
within a lock-in range [0,cl?/) ; the corresponding acquisition 
process in the loop leads, if it is not interrupted, to a new 
locked state without cycle slipping. 

Finally, our definitions give f) loc k-in C fi pu ii-in C ^hold-in- 
If there is a certain stable equilibrium varies continuously 
when cd£ ee is changed within the hold-in, pull-in, and lock- 
in ranges (see Footnote [To|), then 

[0,wz) C [0 ,cj p ) C [0 ,u h ) 

which is in agreement with the classical consideration (see, 
e.g. (67] p.34], (6§ p.612], 0 p.61], [66| p.138], § p.258]). 


D. Approximations of the lock-in range of the classical PLL 

For the case of the classical odd PD characteristic (see 
Fig. 10), taking into account that equilibria are propor¬ 
tional to the frequency deviation (see ( [TT] )) and using the 
symmetry (x eq (ui),0 eq (ut)) = -(x eq (-ui),0 eg (-ui)), we 
can effectively determine cj/. For that, we have to increase 
the frequency deviation |c<;^ ee | step by step and at each 
step, after the loop achieves a locked state, to change 
uj% ee = uj abruptly to cd^ ee = — ui and to check if the loop 
can achieve a new locked state without cycle slipping. If 
so, then the considered value |co^ ee | belongs to Di oc k-in* If 
CJ fr ee _0 belongs to fi pu ii-m, then it is clear that 0 belongs 
to f2i oc k-in (see Fig. [l0| left). The limit value ui is defined 
by the case in Fig. |10[ middle. At the next step when a 
value |c^ ee | = \uj\ > uoi is considered, for cd^ ee = — \uj\ the 
trajectory from the initial point, corresponding to a stable 

right: red trajectory 


free 


= |a;| (see Fig. 


10 


equilibrium for uo A 
outgoing from a black dot), is attracted to an equilibrium 
only after cycle slipping. In other words 77 , for this case: 
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The lock-in range is a subset of the pull-in range such 
that for each corresponding frequency deviation the lock- 
in domain (i.e. a domain of the loop states, where fast 
acquisition without cycle slipping is possible) contains both 
symmetric locked states (i.e. locked states for the positive 


- (red) negative (D A free 

-(black) positive cd a ee 



Fig. 11. Phase portrait. Separatrices, equilibria and corresponding 
local lock-in domains (shaded): upper black is for cn^ ee = 61.5, lower 
red is for = —61, 5. The uniform lock-in domain is approximated 
by the band between two blue horizontal lines: |®| < 0 . 0110 . 


Remark 8. If the filter - perfect integrator can be imple¬ 
mented in considered architecture, the loop can be designed 
with the first order PI filter having the transfer function 
H ( s ) = 1 ~^ T2 . Equations of the loop in this case become 

- (red) negative co A free 


-(black) positive co A ree 



Fig. 12. Phase portraits for the classical PLL with the following pa¬ 
rameters: H(s) = 1 -(°QQ 33 g S 7 L = 250, and o;^ ee = ±47. Separatrices, 
equilibria and corresponding local lock-in domains (shaded): upper 
black is for = 47, lower red is for cj^ ee = —47. The uniform 

lock-in domain is approximated by the band between the two blue 
horizontal lines: |x| < 0.0119. 


and negative value of the difference between the reference 
frequency and the VCO free-running frequency). 

In Fig. [To! middle the set Di ock _ in : contains all equilib¬ 
ria x eq (ujff^) for 0 < |co^ ee | < uji. However for some 
non-equilibrium initial states from the band defined by 
{x : \x\ < \x eq (uji)\} (phase error 6 a takes all possible 
values), cycle slipping can take place. For example, see the 
points to the left and to the right of the black equilibrium 
states (i.e. for cj^ ee = \lji\ > 0), lying above the red 
separatrix (i.e. for cj^ ee = —\wi\ < 0), correspond to the 
red trajectories (i.e. for u;^ ee = — \ui\ < 0), which are 
attracted to an equilibrium only after cycle slipping. To 
approximate the Di ock _ in by a band, ui can be slightly 
decreased to cut the above points. In Fig. [Tl] the band 
defined by Xi ock _ in = {x : \x\ < \x eq (uji)\, < u{\ is 
contained in Di ock _i n and for any initial state from the band 
the corresponding acquisition process in the loop leads, if 
it is not interrupted, to lock up without cycle slipping. 
Such a construction is more laborious and requires rigorous 
analysis of the phase space or exhaustive simulation. 


Remark 7. If we define (see, e.g. [78, p.92]) cycle slipping 
by the interval of maximum length 2i r instead of 47r in 
Definition ^ }* he. limsup^^ |#a(0) — #a(^)| > tt ; then for 
any \cv^ ee \ > 0 a distance between neighboring unstable and 
stable equilibria and a phase deviation of the corresponding 
unstable saddle separatrix may exceed tt (see, e.g. Fig. 11). 
Thus, the lock-in range may contain only |cj^ ee | = 0. 


x = 6 A = ^a 66 - Lx - L—lp(6a ), 

Ti Ti 


(31) 


or equivalently 


e A = -L-ip(e A ) - L-<p\o A )d A . (32) 

Ti Ti 

Here the equilibria are defined from the equations 
f(0eq) = 0, X eq = 

Because model (32) does not depend explicitly on ± ee , 


the hold-in and pull-in ranges are either infinite or empty. 
Note, that the parameter uo^ ee shifts the phase plane ver¬ 
tically (in the variable x) without distorting trajectories, 
which simplifies the analysis of the uniform lock-in domain 
and range (see Fig. [7]| ). If the transfer function H(s) of 
a high order filter has the term s r with r G N in the de¬ 
nominator, then instead of equilibria we have a stationary 

—<jJ ree 

linear manifold: c p(0 eq ) = 0, cix\ q + ... + c r x r eq = —. 

For the classical PLL with the filter’s transfer function 
H( s ) = it can be analytically proved that the pull- 

in range is theoretically infinite. Some needed explanations 
are given by Viterbi W using phase plane analysis. But, 
even in such a simple case, rigorous phase plane analysis 
is a complex task (e.g. [7^, the proof of the nonexistence 
of heteroclinic and first-order cycles is omitted in W- 
The rigorous analytical proof can be effectively achieved 
by considering a special Lyapunov function 


by c 

B 


Jree 


v(x,e A ) = \{ x - u -±-y + 


f sin 2 ^ 


gA > 
2 — 


0 and 
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V(x,0 A ) = — h/3sm 2 6 A < 0. Here it is important that 
for any uj^ ee the set V(x,0 A ) = 0 does not contain the 
whole trajectories of system ( [31] ) except for equilibria. 

E. Initial and free-running frequencies of VCO 

Note that in the above Definitions m and [6] the hold- 
in, pull-in, and lock-in sets are defined by the frequency 
deviation, i.e. by the absolute value of the difference 
between VCO free-running frequency (in the open loop) 
and the input signal’s frequency: \u% ee \ = \ui — ^4 ree |- 
The VCO free-running frequency u4 ree is different from the 
VCO initial frequency c^2 (0): ^(0) = o4 ree + g(0), where 
g(0) = c*x(0) + hip(6&(0)) is the initial control signal, 
depending on the initial states of the filter x(0) and the 
initial phase difference 0 a (0). 

It is interesting that for simplified model (Jl2| with 
h = 0 (see eq. 2.20 in the classic reference [4]) the 
absolute value of the initial difference between frequencies 
|0 A (O)| = |cja( 0)| = |coq — c^2(0) | is equal to the frequency 
deviation || = \uji — cd>2 ree I • Following such simplified 
consideration in engineering literature the concept of an 
“ initial frequency difference ” can be found to be in use 
instead of the concept of a “ frequency deviation see, e.g. 
[3] p.44] “If the initial frequency difference (between VCO 
and input) is within the pull-in range, the VCO frequency 
will slowly change in a direction to reduce the difference ”, 
[80 . p. 1792] u The maximum frequency difference between 
the input and the output that the PLL can lock within one 
single beat note is called the lock-in range of the PLL ”, 
[1, p.49] “ Whether the PLL can get synchronized at all or 
not depends on the initial frequency difference between the 
input signal and the output of the controlled oscillator.” In 
general, the change of cj2 ree to CJ2(0) may lead to wrong 
results in the above definitions of ranges because in the 
case of x(0) 7^ 0, h 7^ 0 or non-odd function p(0 A ) for the 
same values of cj2(0) the loop can achieve synchronization 
or not depending on the filter’s initial state x(0), the initial 
phase difference 6a (0), and Ci4 ree - See the corresponding 
example. 

Example 3. Consider the behavior of model (|8| for 
the sinusoidal signals (i.e. <p(0 A ) = ^sm(26 A )) and the 
fixed parameters: co A = 100, iif(s) = 1 Tl = 

0.0448, t 2 = 0.0185, L = 250. In Fig. [73] the phase portrait 
of system © is shown. The blue dash line consists of 
points for which the initial frequency difference is zero: 
cja(0) = 0a(O) = 0. Despite the fact that the initial 
frequency differences of all trajectories outgoing from the 
blue line are the same (equal to t)), the green trajectory 
tends to a locked state while the magenta trajectory can 
not achieve this. 


VI. CONCLUSIONS 

This survey discussed a disorder and inconsistency in 
the definitions of ranges currently used. An attempt is 
made to discuss and fill some of the gaps identified between 
mathematical control theory, the theory of dynamical sys¬ 
tems and the engineering practice of phase-locked loops. 



Fig. 13. Phase portrait for cj^ ee = 100. Blue dash curve corresponds 
to the set defined by 6a (0) = 0. Initial points of the green (upper) and 
magenta (lower) trajectories correspond to the same initial frequency 
difference u;a(0) «= 0. 


Rigorous mathematical definitions for the hold-in, pull-in, 
and lock-in ranges are suggested. The problem of unique 
lock-in frequency definition, posed by Gardner 19], is solved 
and an effective way to determine the unique lock-in 
frequency is suggested. 
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